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The  Effect  of  Shear  Stresses 


1 


on  the  Carrying  Capacity  of  I-Beams' 

by 

Carl-Fredrik  A.  Leth^ 


Abstract, 

This  paper  considers  the  bending  of  a cantilever  I-beam 
of  a ductile  metal,  such  as  mild  steel,  that  is  loaded  by  a 
transverse  force  at  the  free  end.  Taking  the  length  of  the  beam 
into  account  the  load  carrying  capacity  can  be  expressed  in  terms 
of  the  moment  at  the  built-in  end  that  is  in  equilibrium  with  the 
maximum  transverse  force  at  which  collapse  of  the  beam  is  immi- 
nent, owing  to  the  development  of  regions  of  plastic  flow.  The 
influence  of  the  length  of  the  beam  upon  this  limit  moment  is 
studied  in  the  paper. 

As  the  simplest  approximation  of  the  limit  moment  the 
fully  nlastic  moment  can  be  taken.  From  experimental  evidence 
it  is  known  that  sufficiently  long  beams  can  support  a force 
that  is  in  equilibrium  with  the  fully  plastic  moment.  But  it  is 
also  observed  that  for  beams  that  have  both  a short  length  and 
an  I-shaped  cross-section,  the  limit  moment  may  be  considerably 
less  than  the  fully  plastic  moment.  The  aim  of  the  paper  is  to 
obtain  theoretical  estimates  of  the  maximum  moment  that  these 
I-beams  can  develop  at  the  built-in  end,  assuming  that  the 
material  exhibits  ideal  plasticity, 

1.  The  results  presented  in  this  paper  were  obtained  in  the 
course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  N7onr-35^01  with  Brown  University, 

2.  Research  Assistant,  Graduate  Division  of  Applied  Mathematics, 
Brown  University,  Providence  12,  R,  I, 
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At  present  a complete  three-dimensional  solution  of  the 
problem  seems  to  be  out  of  the  question.  Even  for  a two-dimen- 
sional model  of  the  problem  it  would  be  very  difficult  to  obtain 
a complete  elastic-plastic  solution.  The  limit  analysis  theorems 
of  Drucker,  Prager  and  Greenberg  are  therefore  applied.  These 
furnish  upper  and  lower  bounds  for  the  limit  moment.  Two  approx- 
imate lower  bounds  are  obtained  by  constructing  two  different 
admissible  stress  fields  for  a two-dimensional  model  of  the 
problem.  Only  the  second  bound  gives  estimates  of  the  limit 
moment  for  short  beams,  but  both  bounds  are  applicable  to  long 
beams.  For  short  beams  an  upper  bound  is  obtained  from  a kine- 
matically admissible  velocity  field,  using  the  appropriate  limit 
analysis  theorem. 

The  first  stress  field  is  based  on  a natural  extension 
of  the  conventional  elastic  theory  of  the  bending  of  beams, 
modifying  this  in  such  a way  that  at  loads  above  the  load  at 
which  plastic  flow  first  occurs  the  yield  condition  is  not 
violated  anywhere.  The  second  stress  field  is  artificial  in  the 
sense  that  it  is  unrelated  to  the  stress  distribution  of  an 
elastic  beam.  Nevertheless,  it  furnishes  in  a much  shorter  way 
than  the  first  a lower  bound  for  the  limit  moment,  valid  for 
both  short  and  long  beams. 

Both  stress  fields  yield  values  of  the  limit  moment  which 
tend  towards  the  fully  plastic  moment  for  long  beams  as  the  ratio 
of  the  beam  length  to  depth  is  increased.  For  short  beams  the 
second  stress  field  gives  a value  of'  the  limit  moment  which 
differs  considerably  from  the  fully  plastic  moment.  For  this 
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reason  a velocity  field  is  constructed  so  that  an  upper  bound 
is  obtained  for  short  beams.  This  velocity  field  is  based  on 
the  failure  of  the  short  I-beam  primarily  because  of  shear  in 
the  v/eb.  It  is  shown  by  examples  that  the  upper  and  lower 
bounds  are  very  close  for  short  beams# 

The  theoretical  results  obtained  in  this  paper  are  com- 
pared with  the  reported  experimental  data  that  are  concerned 
with  the  influence  of  shear  stresses#  The  experiments  give  a 
reasonable  confirmation  of  the  theoretical  results  both  for  the 
value  of  the  limit  moment  and  for  the  computed  critical  length 
of  the  beam  below  which  the  carrying  capacity  rapidly  decreases, 

A.  Introduction. 

The  problem  of  the  elastic-plastic  bending  of  beams  is 
so  complex  that  at  present  no  rigorous  complete  solutions  have 
been  found.  In  this  paper  is  presented  an  approximate  treat- 
ment of  a cantilever  I-beam  which  carries  a transverse  force 
parallel  to  the  web  at  the  "free"  end,  Fig,  1,  The  purpose  is 
to  study  the  effects  of  the  shear  stresses  on  the  carrying 
capacity  of  the  beam.  These  may  gain  importance  because  of  the 
difference  in  width  of  the  flange  and  the  web. 

In  the  elastic  analysis  of  the  general  problem  of  a 
beam  under  end  loading  an  exact  solution  can  be  obtained  by  the 
seim-inverse  method  of  Saint-Venant , For  commonly  used  beam 
sections,  however,  this  solution  is  too  complicated  to  be  used 
in  normal  engineering  practice.  Instead  an  approximate  solution 
is  used  that  is  based  upon  an  assumption,  introduced  by  J« 
Bernoulli,  concerning  the  deformation  of  the  beam.  The  deflection 
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of  the  beam  is  assumed  to  take  place  in  such  a way  that  points 
which  initially  lie  in  a plane  normal  to  the  center-line  after 
the  deformation  also  form  a plane  that  is  normal  to  the  de- 
flected center-line.  This  condition  together  with  the  equili- 
brium equations  and  Hooke1 s law  yields  results  that  do  not  ful- 
fill the  compatibility  equations.  Furthermore  if  the  cross- 
section  is  not  rectangular  the  condition  of  a stress-free 
cylindrical  surface  of  the  beam  is  violated.  These  inconsis- 
tencies are  generally  of  minor  importance,  however,  as  many 
experiments  have  snown. 

In  considering  the  bending  of  the  beam  of  Fig.  1 beyond 
the  elastic  limit,  we  have  formulated  a problem  that  is  of 
fundamental  importance  in  the  theory  of  limit  analysis  of  struc- 
tures. It  is  not  appropriate  in  this  paper  to  discuss  the 
procedures  of  limit  analysis.  Such  discussions  can  be  found 
e.g.  in  the  book  by  Van  den  Broek  [1]  and  the  papers  by  Baker 
[2]  and  Symonds  and  Neal  [3],  It  is  sufficient  to  state  that 
limit  analysis  of  beams  and  frames  is  an  analysis  of  failure  of 
these  structures  associated  with  the  formation  of  "plastic 
hinges"  in  the  members  of  the  structure.  Plastic  hinge  sections 
are  defined  as  those  across  which  rotations  of  arbitrary  magni- 
tudes can  occur  while  a constant  bending  moment,  called  the 
limit  moment,  M0,  is  transmitted.  In  this  paper  we  discuss  the 
computation  of  the  hinge  moment  at  the  base  of  a cantilever  I- 
beam,  since  this  can  be  regarded  as  the  basic  structural  element 

* Numbers  in  square  brackets'  refer  to  the  Bibliography  at  the 
end  of  the  paper. 
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in  any  continuous  beam  or  rigid  frame.  We  shall  compare  the 
limit  moment  for  such  a cantilever  with  the  approximation  to  the 
limit  moment  that  is  commonly  used  in  the  applications  of  limit 
analysis,  namely  the  "fully  plastic  moment",  M . This  latter 
moment  is  now  considered  in  detail. 

Consider  first  a beam  subjected  to  pure  bending.  Let 
the  beam  have  two  planes  of  symmetry  and  let  the  moment  act  in 
one  of  them,  Fig,  2,  Assume  the  material  to  follow  the  perfect- 
ly elastic-plastic  stress-s train  law,  i.e,,  the  stress-strain 
diagram  in  simple  tension  of  Fig.  3*  Moreover  let  it  be  assumed 
that  each  fiber  is  in  a state  of  simple  tension  or  compression, 
so  that  d is  the  only  non-vanishing  stress.  The  condition  that 
cross-sections  remain  plane  in  the  elastic  part  of  the  beam  then 
yields  the  distribution  of  stress  over  the  cross-section  that  is 
depicted  in  Fig.  In  Fig,  4a  the  outmost  fiber  has  not  reach- 

ed the  yield  stress,  but  in  4b  the  moment  has  increased  so  that 
the  outer  part  of  the  cross-section  becomes  plastic,  while  the 
inner  part  remains  in  the  elastic  state.  In  Fig,  4c  the  curva- 
ture of  the  beam  has  taken  such  a large  value  that  the  elastic 
core  in  the  center  can  be  neglected  and  the  plastic  regions  are 
considered  to  cover  the  entire  cross-section.  The  bending  moment 
corresponding  to  this  last  state  is  called  the  fully  plastic 
moment • 

In  general  for  a beam  having  two  planes  of  symmetry,  in 
one  of  which  the  load  is  applied,  the  fully  plastic  moment  is 
given  by 

mp  = YZP 
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where  Z is  twice  the  moment  of  the  cross-section  area  on  one 
side  of  the  middle  surface,  with  respect  to  the  transverse  axis 
of  symmetry.  In  the  special  case  of  an  I-beam,  the  fully  plastic 
moment  can  be  considered  to  consist  of  two  parts,  namely  the 
contributions  from  the  flanges  and  from  the  web.  We  then  have 


+ M 

pw 


where  M = b, (a2  - c2)Y  and  M = bc2Y. if  the  web  and  flange 
pf  1 pw 

are  assumed  to  be  rectangles,  Fig.  5. 

Note  that  the  stress  distributions  of  Fig.  4 are  not  an 
exact  solution  even  for  a rectangular  beam  in  pure  bending,  as 
was  pointed  out  by  Hill  [4];  these  stresses  are  not  consistent 
with  continuity  of  displacements  at  the  elastic-plastic  boundary. 

In  the  more  general  problem  of  the  bending  of  the  beam 
in  the  presence  of  shear  forces,  the  moment  varies  along  the 
beam,  and  it  is  necessary  to  remove  the  previous  assumption  that 
t is  zero.  It  is  shown  in  the  book  of  Prager  and  Hodge  [5]  , 
see  also  Hodge  [6],  that  must  vanish  in  the  plastic  region 

that  spreads  in  from  the  outer  fibers.  Hence  the  distribution 
of  the  longitudinal  stress,  depicted  in  Fig.  4,  remains  unchanged 
if  the  other  assumptions  of  the  analysis  for  simple  bending  are 
retained.  The  depth  of  the  plastic  region  then  varies  according 
to  the  distribution  of  the  moment  along  the  beam.  Figure  6 shows 
the  plastic  region  in  the  cantilever  beam  loaded  with  a shear 
force  at  the  free  end.  As  the  elastic  region  must  carry  the 
total  shear  force,  the  shear  stresses  increase  in  magnitude 
towards  the  built-in  end  at  the  center  of  the  beam.  Failure  by 
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collapse  occurs  when  the  limit  moment  in  the  presence  of  shear 
is  reached  at  the  built-in  sections;  deflections  would  then 
continue  under  constant  load,  if  geometry  changes  are  neglected, 
with  rotations  of  arbitrary  amounts  at  the  plastic  hinge  at  the 
built-in  end.  Tn  most  applications  of  limit  analysis  the  com- 
ponent beams  are  assumed  to  be  sufficiently  long  with  respect 
to  their  depth  so  that  the  shear  stresses  can  be  neglected  in 
computing  the  limit  moment.  Thus  the  fully  plastic  moment  is 
taken  in  these  cases  as  a good  approximation  of  the  actual  limit 
moment . 

An  early  paper  dealing  with  the  shear  stresses  at  the 
section  of  failure  of  the  beam  was  written  by  Stu’ssi  [7].  He 
considers  a beam  under  the  influence  of  a moment  and  a shear 
force.  Assuming  that  the  longitudinal  stress  follows  a general 
stress-strain  law  in  simple  tension  for  mild  steel  and  consider- 
ing plane  cross-sections  to  remain  plane  over  the  entire  cross- 
section,  the  distribution  of  shear  stress  is  obtained  with  the 
help  of  the  equilibrium  equation.  The  fact  that  the  value  of 
the  longitudinal  stress  at  yield  is  affected  by  the  presence  of 
these  shear  stresses  is  neglected.  The  shear  stress  at  the 
center  of  the  critical  cross-section  at  a specified  stage  of  the 
loading  is  considered  to  form  the  condition  for  the  failure  of 
beam. 

A recent  investigation  of  the  shear  stresses  at  the  most 
highly  stressed  section  has  been  presented  by  Horne  [8],  A 
rectangular  or  an  I-shaped  cross-section  under  the  assumption 
of  plane  stress  and  the  Tresca  yield  condition  is  treated.  At 
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failure  two  different  plastic  regions  appear.  One  of  these 
contains  the  outer'  fibers  of  the  beam  at  the  built-in  end  where 
longitudinal  stresses  are  large 5 the  other  plastic  region  devel- 
ops at  the  center  of  the'  beam  where  shear  stresses  predominate. 
The  assumption  that  plane  cross-sections  remain  plane  during  the 
deformation  is  adopted  in  the  elastic  region.  In  the  plastic 
regions  account  is  taken  of  the  equilibrium  and  yield  conditions, 
but  not  of  strains.  It  is  not  possible  in  this  solution  to 
satisfy  all  the  matching  conditions  of  the  stresses  at  the 
plastic-elastic  boundaries.  For  the  case  of  I-beams  the  poss-. 
ibility  of  a violation  of  the  yield  condition  at  the  flange-web 
junction  is  not  investigated. 

A different  viewpoint  is  presented  in  a paper  by  Onat 
and  Shield  [9].  The  same  problem  as  in  Fig.  1 is  treated  for  a 
rectangular  cross-section  with  the  assumptions  of  plastic-rigid 
material  and  plane  strain  conditions.  The  exact  solution  (ful- 
filling both  stress-  and  flow-conditions)  for  the  region  near 
the  built-in  end  is  obtained. 

In  the  experimental  studies  of  the  bending  of  beams 
beyond  the  elastic  limit,  there  are  many  results  showing  that 
for  long  -beams  the  fully  plastic  moment  provides  a good  approxi- 
mation to  the  limit  moment.  A survey  of  these  studies  has  been 
presented  by  Roderick  and  Phillips  [10].  On  the  other  hand, 
relatively  few  experiments  have  been  made  whose  aim  was  to  study 
the  behavior  of  the  beam  when  shear  effects  are  important.  - 
Apparently  the  only  extensive  tests  of  this  nature  are  those 
made  by  Baker  and  Roderick  [11],  Hendry  [12],  and  Johnston, 
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Yang  and  Beedle  [13].  The  comparison  of  the  results  of  this 
paper  -with  these  experiments  is  postponed  until  after  the  analy- 
sis has  been  described. 

B.  General  concents. 

The  analysis  in  this-  paper  is  based  upon  the  two  limit 
analysis  theorems  of  Drucker,  Prager  and  Greenberg  [ l1!]  » The 
main  part  of  the  paper  is  concerned  with  the  first  of  these 
theorems,  and  we  begin  by  defining  the  terms  and  concepts  that 
are  used  in  it, 

A yield  function  3>  is  a positive  definite  function  of 
the  stress  components  only,  such  that  plastic  flow  can  occur 
when  <5  = 1.  In  a stress-free  state  we  have  = 0 and  the  elas- 
tic range  consists  of  0 < $ <1.  We  then  define  a safe  state  Of 
stress  as  a state  for  which  # < 1 throughout  the  structure. 

On  the  other  hand  collapse  is.  defined  as  the  state  for  which 
plastic  flow  would  occur  under  constant  load  if  the  accompanying 
changes  in  the  geometry  of  the  structure  were  disregarded..  More- 
over we  call  a stress  system  statically  admissible  when  it 
satisfies  the  equilibrium  equations  and  is  consistent  with  the 
surface  tractions  on  the  boundary  of  the  body»  The  first  limit 
analysis  theorem  then  states  that  if  a safe  statically  admissible 
state  of  stress  can  be  found  at  each  stage  of  loading,  collapse 
will  not  occur  under  the  given  loading  schedule . Thus  a load 
at  which  a safe  statically  admissible  stress  field  can  be  con- 
structed is  a lower  bound  for  the  collapse  load  of  a structure. 

In  the  main  part  of  this  paper  we  apply  the  theorem 
stated  above  to  the  problem  of  estimating  a lower  bound  for  the 
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maximum  end  force  that  can  be  applied  to  a cantilever  beam,  Fig. 
1.  As  discussed  later  in  this  section  we  find  it  necessary  to 
make  certain  approximations  in  the  construction  of  safe  stati- 
cally admissible  stress  fields  fur  this  problem.  Since  the 
conditions  of  the  theorem  are  not  satisfied  exactly,  the  maximum 
loads  obtained  in  our  analysis  cannot  be  regarded  as  true  lower 
bounds  of  the  collapse  load  of  the  problem.  However,  it  is 
believed  that  the  nature  of  the  approximations  is  such  that  the 
results  will  be  of  practical  usefulness,  and  will  serve  as  a 
guide  to  further  experiments. 

As  it  is  part  of  our  purpose  to  study  plastic  regions  in 
the  beam  we  consider  stress  systems  that  fulfill  the  yield  con- 
dition $ < 1,  the  equilibrium  equations,  and  agree  with  the 

specified  surface  tractions.  We  then  know  that  the  correspond- 
ing load  either  is  the  true  collapse  load  or  is  smaller  than 
this  load.  This  is  true  because  of  the  fact  that  any  smaller 
load  will  be  in  equilibrium  with  a safe  statically  admissible 
stress-system  and  therefore  cannot  be  a collapse  load.  Note 
that,  when  using  this  theorem,  we  focus  our  attention  entirely 
on  the  stresses.  This  was  done  also  by  Horne  [8],  who  dealt 
primarily  with  conditions  on  the  stresses.  The  lower  bound 
theorem  provides  a basic  justification  for  this  point  of  view. 

We  shall  here  use  the  yield  function  that  was  proposed 
by  von  Miscs.  It  has  the  form  (see  [5]) 


+ °y2  + t2 


d d - d cJ  -q  q 
x y y z z x 
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Although  we  have  eliminated  above  all  considerations  of 

compatibility  of  strains  the  problem  is  still  so  complicated  that 

further  simplifications  must  be  made.  Following  the  method  of 

the  approximate  elastic  solution,  mentioned  in  the  Introduction, 

this  analysis  is  made  two-dimensional  in  the  x,z  plane,  so  that 

d , t and  d^  are  the  only  non-vanishing  stresses.  Moreover 
x xz  z 

we  assume  the  stresses  to  be  independent  of  the  y-coordinate . 

This  implies  that  the  condition  of  no  surface  traction  on  the 
inner  sides  of  the  flanges  must  be  violated.  The  same  error 
appears  in  the  approximate  elastic  solution.  This  cannot,  of 
course,  be  taken  as  a justification  of  the  assumption  in  the 
problem  of  the  plastic  bending  of  the  beam.  Such  justification 
must  ultimately  be  found  in  a more  complete  theory  or  by  experi- 
ments. 

Two  different  stress  fields  are  constructed  in  the 

following.  In  the  first  analysis  (Section  C)  it  is  desirable  to 

make  a further  approximation.  This  consists  of  assuming  the  d„ 

z 

stress  to  vanish,  and  it  turns  out  that  one  of  the  two  equili- 
brium equations  cannot  generally  be  satisfied  in  this  case.  In 
the  second  analysis  (Section  D)  the  same  assumption  is  made,  but 
for  this  stress  field  it  is  found  that  a vanishing  d^  docs 
satisfy  both  equilibrium  equations.  Hence  no  approximation  in 
this  regard  is  involved  in  the  second  type  of  stress  field.  As 
this  second  analysis  fives  quite  satisfactory  values  of  the 
estimates  of  the  limit  moment  no  attempt  is  made  to  consider  the 

effects  of  the  a stress  in  the  first  analysis. 

z 
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The  first  analysis  is  built  up  as  a continuation  of  the 
conventional  elastic  theory  of  the  bending  of  beams.  In  this 
theory  it  is  known  that  the  most  highly  stressed  points  of  long 
beams  are  located  in  the  outer  fibers  at  the  built-in  section. 
Hence  we  assume  that  plastic  regions  start  at  these  points  and 
then  spread  in  through  the  flange  into  the  web,  Fig.  6,  At  some 
stage  a new  plastic  region  will  form  starting  from  the  center 
point  of  the  web  at  the  built-in  section.  In  this  first  analy- 
sis we  estimate  the  limit  load  by  the  load  that  corresponds  to 
the  stage  ween  the  stresses  at  the  center  of  the  beam  at  the 
built-in  section  just  start  to  produce  plastic  flow.  The  present 
estimates  of  limit  loads  are  conservative  on  this  account.  The 
calculation  in  the  first  part  of  the  analysis  considers  this 
case  where  the  plastic  regions  have  the  shape  indicated  in  Fig, 

6.  However,  when  the  yield  condition  is  checked  in  the  elastic 
region,  we  find  that  a critical  point  exists  at  the  flange-web 
junction  for  beams  of  relatively  short  length.  This  suggests 
that  a second  plastic  region  will  occur  in  the  outer  portions 
of  the  web,  as  shown  in  Fig,  7»  The  second  part  of  this  analysis 
considers  this  case.  In  the  third  part,  beams  of  still  shorter 
lengths  are  treated,  namely  the  case  when  the  second  plastic 
region  extends  into  a cross-section  of  the  beam  that  otherwise 
would  have  been  entirely  elastic.  This  is  shown  in  Fig.  8. 

At  the  start  of  the  second  analysis  (Section  D)  a stress 
field  is  constructed  which  is  anticipated  by  the  failure  of  a 
very  short  beam  due  to  shear  in  the  web.  It  is  then  found  possi- 
ble to  generalize  this  stress  field  so  that  beams  of  arbitrary 
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lengths  can  be  treated. 

The  second  type  of  stress  field  is  both  simpler  and  more 
general  than  the  first  type,  and  it  will  probably  provide  the 
more  useful  basis  for  estimating  limit  moments  in  the  presence 
of  shear  forces.  On  the  other  hand  the  second  type  of  stress 
field  is  quite  artificial,  whereas  the  first  type  is  a natural 
extension  of  the  veil  confirmed  simple  theory  of  elastic-plastic 
bending  of  beams.  In  this  respect  it  corresponds  to  and  supple- 
ments the  analysis  of  Horne  [8] . 

C.  Analysis  for  the  first  type  of  stress  field. 

For  the  first  type  of  stress  field  this  analysis  is 
divided  into  three  parts,  based  in  turn  on  the  three  configura- 
tions of  plastic  regions  shown  in  Figs.  6-8,  The  analysis  will 
show  that  there  exist  for  a given  cross-section  shape  three 
decisive  lengths  £ and.  which  determine  the  appropriate 
configuration. 

1,  Beams  of  length  so  that  I > £ ^ and  ^ 

Vie  here  consider  plastic  regions  as  shown  in  Fig.  6. 

They  are  bounded  by  the  outer  sides  of  the  flanges,  by  the  built- 
in  section  and  the  curve  u(:c).  In  zone  BC  of  the  beam  the 
plastic  regions  extend  into  the  web,  but  in  zone  CF  the  bounda- 
ries of  the  plastic  regions  are  located  in  the  flange.  With  the 
assumption  = 0,  the  two  remaining  stresses  in  the  plastic 
regions  are  described  by  the  equilibrium  equation  and  by  the 
yield  condition: 
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X 
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xz 
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(2) 

The  yield  condition  (2)  is  identically  satisfied  if  we  relate 
the  stresses  to  each  other  through  the  function  a(x,z)  in  the 
following  manner. 


dx  = Ycosa,  Txz=-|Sinai 

The  equilibrium  equation  becomes 

-|3  sin  a - cos  a = 0.  (3) 

ox  az 

The  characteristics  of  this  equation  are  lines  along  which  a is 
constant.  The  expression 

da  = H + ^ = ° 

together  with  equation  (3)  implies  that  the  characteristics  are 
straight  lines  with  the  slope 

cot  a.  (4) 

On  the  boundary  z = a we  have  d = Y,  % =0  or  a(x,a)  = 0. 

x 7 xz  1 

The  characteristics  are  therefore  lines  parallel  to  the  z-axis 
and  the  stresses  in  the  plastic  region  are  given  by 

cx  = Y>  Txz  = °* 

In  the  elastic  part  of  the  beam  we  follow  the  theory 
based  upon  the  Bernoulli  assumption.  This  implies  a tensile 
stress  linear  in  z.  In  the  elastic  part  of  zone  BC  we  therefore 
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have 


d 


x 


Y. 


(From  here  on  we  consider  only  positive  z-values.  The  quantities 

in  the  beam  are  either  symmetrical  or  anti-symmetrical  around  the 

midale  surface.)  The  equilibrium  equation  (1)  determines  the 

shear  stress.  Using  the  boundary  condition  t (x,u)  = 0 we 

obtain  , 2 

^xz  = - \ (1  “ ^)Y 

u 

where  the  prime  denotes  differentiation  with  respect  to  x*  The 
boundary  u(x)  between  the  plastic  and  elastic  regions  in  the  zone 
BC  is  determined  by  the  total  moment  at  any  cross-section.  It 


is 


or 


Px  = mp  - f “s  M 


pw 


u2  M - Px 

S = 3 . 

c pw 

The  cross-section  C is  determined  bv  uc  = c or 


XC  = 


M - i M 
P 3 pw 


(5) 


(6) 


In  the  elastic  part  of  zone  CF  we  have 


d 


x 


Y 


and  the  shear  stress  in  the  flange  is  also  given  by  equation  (1) 

Txz  = \ (1  - Y. 

u 

To  find  the  shear  stress  in  the  web  we  must  consider  the  differ- 
ence in  the  width  of  the  beam  at  the  flange-web  junction.  We 
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let  the  equilibrium  condition  be  satisfied  only  approximately; 
namely  v/e  assume  that  the  shear  stress  has  the  same  jump  at  the 
flange-web  junction  as  the  width  of  the  beam.  Hence  the  boundary 
condition  for  the  shear  stress  in  the  web  is 


' xz 


(x,c)  = - ^ (1  - -4)  Y. 


T- 

u 


The  equilibrium  equation  (1)  then  gives  the  shear  stress  in  the 
web 


' xz 


,2  b-,  2 

= ‘ \ [1  - h + <■ b - na  - %)]y. 

u u 


As  before  the  plastic-elastic  boundary  u(x)  in  zone  CF  is  deter- 
mined by  the  total  moment  at  any  cross-section.  It  is 


Px  = *|  (a2  - u2)bF  + (u2  - c2) 


£ b, 
u l 


+ (u.  c).2  bi[0  + 2(u  . c)]  + | £ b|Y 


u 


or 


U3  ,(£. 

c3  3 


Pxb  i u , . b_ 


M^,,b, 


3S  + 2d  - = o. 


■pw  1 M1 

The  cross-section  F is  determined  by  u^,  = a,  i.e,, 


(7) 
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bl 


= 0 


or 


x = II 

F aP 


(8) 


where  I is  the  moment  of  inertia  oi  the  cross-section 


= | [b-^3  - (b1  - b)c3] , 


In  the  entirely  elastic  zone  FH  the  tensile  stress  is 


given  by 
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3X  = f XZ. 


The  equilibrium  equation  (1)  again  determines  the  shear  stress, 
so  that 

P(a2  - z2) 


xz 


21 


in  the  flange 


M 


C2  - Z2) 

M 

t - ^ in  the  web. 

xz  2I 

The  validity  of  the  stress-system,  obtained  above, 
depends  on  two  conditions.  First,  it  is  assumed  in  ^ig.  6 that 
the  plastic  regions  extend  into  the  web  at  the  built-in  section, 
Secondly,  the  yield  function  is  assumed  not  to  exceed  unity  in 
the  elastic  region.  The  first  condition  implies  Ug  c,  i.e., 


I > 


M - r M 
P 3 pw 


(9) 


For  the  second  condition  we  form  the  yield  function 

® = (o!  + 3tx!)  ^ 

in  the  different  parts  of  the  elastic  region.  We  obtain 

-2  on  n*  j—  2 


4> 


BC 


= + 27  /£c_v2  c2  n z2  v2 

..2  + 16  '■M  ' T " T) 

U DW  ,1 


-pw  u 
CF  2 + 16  { 


u 


2 on  n.  2 (u/c)2  - (z/c)2  .2 


pw 


in  the  flange 

[(u/c)3  - l]b]_/b  + 1 

2 1 - (z/c)2  + [(u/c)2  - l]bi/b  2 

<X>  = £ + 27  (££_,2  ( 1 ' 

CF  u2  16  Mpw  [(u/c)3  - l]bx/b  + 1 


in  the  web 
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FH 


<X> 


FH 


[x2z2  + jl  (a2  - z2)^] 


i2y2 


2 o v2 

7U.  1 


^ Cx2z2  + 1 $ 


l2Y2 


- z2)2] 


pw 


in  the  flange 
in  the  web. 


When  we  look  for  possible  maxima  of  <1*  we  notice  that  in  the  zone 
FH  <5  has  its  maximum  value  at  the  section  F.  Hence  we  need  only 
to  check  the  zones  BC  and  CF. 

2 2 

In  the  zone  BC  it  is  convenient  to  consider  z /u  = 5 
and  27/16  (Pc/MpW)2  c2/u2  = q as  the  two  independent  variables. 
Then 

*BCU,n)  = 5 +n  a - O2. 


A maximum  appears  on  the  boundaries  t)  = Tmax  and  5 = 1,  At 
5 = 1 we  have  <5  = 1,  so  the  critical  section  is  r]  = corres- 

ponding to  the  smallest  possible  value  of  u,  i.e.  , the  built-in 
section.  The  maximum  of  <5  at  this  section  appears  for  z = 0. 

As  mentioned  in  the  previous  section  we  take  the  load  as  an 
estimate  of  the  limit  load  in  this  analysis,  where  P]_  is  deter- 
mined by  the  yield  condition  at  this  critical  point.  Hence  we 
have 


on  Pi  C n 

ft  <rL'2 


'M  ’> 
pw 


= 1 


u 


B 


u is  eliminated  by  means  of  equation  (5). 
B 


This  gives 


8 %?_  MPW  r 

9 C2  MP 


1 + 


2 M 

C -2- 

7^  M 

b pw 


- 1]. 


(10) 


The  yield  condition  remains  to  be  checked  in  zone  CF. 
We  find  when  looking  for  extreme  values  that  there  are  no 
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stationary  points  in  the  flange.  On  the  boundaries  of  the  elas- 
tic region  in  the  flange  we  have  on  one  side  the  plastic  region 
and  on  the  other  side  the  junction  to  the  web.  Because  of  the 
jump  of  the  shear  stresses  at  this  junction  the  stresses  in  the 
web  are  more  critical.  The  yield  function  in  the  web  also  has 
maximum  on  the  boundaries,  i.e.,  z = 0 and  z = c.  At  z = 0 we 
have 

,\2 


$ 


CF 


- a 


16 


pw 


[(u/c)*"  - 1 ]b^/b  + IVs 
[(u/c)3  - l]b1/b  + 1 


As  u/c  > 1 $>  has  maximum  for  smallest  possible  u/c,  i.e., 

u = c.  But  for  this  value  we  know  from  zone  BC  that  the  yield 
condition  is  satisfied.  Hence  we  have  to  check  at  the  flange- 
web  junction.  We  have  at  z = c 


$ - c 


CF 


2 


u 


/CCu/cT  - l]b  /b  \ 

2 2 (ESL.)2f  l \ 

16  M J I •>  i 

Pw  \^(u/cy  - l]b1/b  + 1 J 


Equating  <3>  = 1 yields  the  following  equation  in  u/c: 

Or 


or 


«3 


r f~2 

- i + JL  = £s_  u /u  . i 

b,  4^  M c \/ 

1 pw  V c 


(11) 


4 


_ 27  _ 2(1  - j2L)(uy 

° ^ SMpw  c Vc 


( 


(H)2  + (i  - ^)2  = 0. 
M 


pw  1 

The  last  equation  shows  that  there  are  two  or  zero  positive 
roots.  We  can  get  an  approximate  solution  by  considering  the 
quantity  u/c  - 1 = e,  which  is  assumed  to  be  a small  number  as 
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compared  to  unity.  The  equation  becomes 


PC3-#  Pc  t -\  , „ >>  . /-  — o . r~  i o ^ _ 2 b 


^ (1  + e ) \/e"  V 2 - \/7  (3  + 3e  + e ) c 

pw  1 


As  ^ 3 \/3  Pc/M  may  be  of  the  same  order  as  \fz  , we  get  in 
approximate  form 

Vl  ( 2-^S  - 3 ✓!)  = r^- 


pw 


or 


e = 


- 1 19  fPc  i2  b , ,(■  9 ,Po  ,2  b i2  . b . 29 

3 Mpw  ' K ~ P Mw)  ' bl  " V j • ^ 


The  critical  value,  when,  yield  just  starts  at  one  point,  corres 
ponds  to  vanishing  of  the  radical,  i.e., 


4 

3 


V 


)2b 

V 


located  at  the  section  where 


(13) 


In  the  cases  of  I-scctions  of  particular  interest  here  b/b^  0.1 
so  that  the  assumption  that  e is  small  as  compared  to  unity  is 
satisfied. 

We  now  consider  t Ac  beam  loaded  with  the  force  P-^  accord 
ing  to  equation  (10),  Then  elimination  of  between  (10)  and 
(13)  gives  an  approximate  value  of  the  length  of  the  beam, 
for  which  the  two  critical  points  start  to  yield  simultaneously. 


(14) 
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On  the  other  hand  eliminating  P1  between  the  inequality  (9), 
taken  as  an  equation,  and  equation  (10)  gives 


pw 


(15) 


Hence  we  have  the  following  two  conditions  for  the  validity  of 
this  first  part  of  the  analysis: 


l 


ir  > and  - > -2  . 

c — c c --  c 


It  is  interesting  to  compare  these  two  values  for  standard  I- 
beams.  Table  1 shows  that  for  many  common  sections  ^ > -lQ* 
This  shows  that  the  stresses  at  the  flange-web  junction  play  a 
decisive  role  even  for  relatively  long  beams. 

The  influence  of  the  shear  stresses  upon  the  limit 
moment  is  expressed  by  equation  (10).  This  influence  is  shown 
in  Fig.  9j  where  P-,£/M  is  plotted  as  function  of  — ,/M  /M  • 

’ 1 P c\/PwP 


2.  Beams  of  length  so  that  _>  Z > and  ^ > -t 0 * 

We  now  consider  beams  that  have  a length  shorter  than 
- 1 ■]_.  We  retain  the  previous  estimate  of  the  limit  load  as  the 
force  according  to  equation  (10).  The  stress  field  must 
therefore  be  modified  so  that  the  yield  condition  $ 1 is 

satisfied  everywhere.  We  introduce  second  plastic  regions  in 
the  web  as  shown  in  Fig.  7»  These  regions  are  bounded  by  the 
curve  v(x)  and  occupy  the  zone  DE,  where  we  assume,  in  this 
part  of  the  analysis,  that  the  cross-sections  D and  E are 
located  between  the  cross-sections  C and  F. 
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There  is  no  reason  to  change  the  previous  stresses  in 
the  zones  BD  and  EH  of  the  beam.  The  cross-sections  D and  E are 
determined  by  equation  (11)  or  approximately  by  equation  (12). 
The  condition  that  the  cross-section  E will  be  located  to  the 
left  of  cross-section  F is  given  by  u^  <_  a.  By  means  of  equa- 
tions (8)  and  (11)  this  condition  is 


where 


K 


2a 

IY 


M M -i 
pw  pf 


(15) 


We  also  assume,  as  before,  that  u^  < c.  The  bounds  for  the 
length  of  the  beam  for  which  this  part  of  the  analysis  applies 
therefore,  are 


and 


5 


Now  consider  zone  DE.  Both  the  flange  and  the  web  are 

partly  elastic,  partly  plastic.  In  the  plastic  region  of  the 

flange  we  have  as  before  d = Y,  t = 0,  In  the  plastic  region 

of  the  web  we  again  relate  the  stresses  to  each  other  by 

dx  = Y cos  a(x,z)  and  %xz  - -X  sin  a(x,z).  If  the  boundary 

/3 

condition  at  the  flange-web  junction  were  given,  we  then  could 
determine  the  stresses  along  the  characteristics.  Therefore  we 
consider  a(x,c)  = f(x)  as  a given  function  and  express  further 
quantities  in  terms  of  it„ 

In  the  elastic  region  of  the  flange  the  assumption  of 
linear  variation  of  dx  ana  of  continuity  of  dx  and  txz  across 
the  elastic  plastic  interface  then  yields 
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6 = (u  - z)  cos  f + z - c y, 

°x  u - c 


The  equilibrium  equation  (1)  together  with  the  boundary  condition 
txz(x,u)  = 0 gives 


xz 


= I1- 


The  other  boundary  condition 

t (x,c)  = — sin  f 

xz  v/3  bx 

determines  the  function  u(x)  to  within  an  arbitrary  constant  A, 
as  follows : 


x^ 
i D 


(u  - c)(l  - cos  f)  = A + 


2b 


* h Jx 


sin  f dx. 


The  boundary  condition  u = uD  at  cross-section  D determines  A 

so  that  I'  ° 

2b  | sin  f dx 

— . (16) 


(uD  - c)(l  - cos  f (::D) ) 
u = c + — + 


1 - cos  f , . . 

s/3  b^Cx  - cos  f) 

But  u is  given  also  at  section  E.  Hence  we  must  impose  one 
condition  on  the  function  f,  namely 

X_ 

v/3  b 2d 

1 [(u^-c) (1-cos  f(xE)  )-(u^-c)  (1-cos  f (xD) ) ] = j sin  f dx. 

Jx^  (17) 

This  equation  can  also  be  derived  as  the  equilibrium  equation 
of  forces  in  the  x direction  on  the  segment  at  the  flange  between 
D and  E. 

In  the  second  plastic  region  in  the  web  the  stresses  are 
given  along  each  characteristic.  Hence  d = Y cos  f(x)  and 

A 


All-10? 


2^ 


= -4o  sin  f(x)  where  the  value  of  x refers  to  the  coordinate 
xz  ^3 

of  the  characteristic  at  the  flange-web  junction,  Fig.  10.  As 
before,  the  characteristics  are  straight  lines  in  this  case 
making  the  angle  <p(x)  to  the  z-direction.  In  view  of  equation 

- 4^  = </3  tan  f - tan  cp*  ' (18) 

Q. 

It  is  found  convenient  to  denote  the  length  of  the  characteris- 
tic by  g(x)  \/l  + 2 sin2f(x). 

Now  consider  a "cross-section"  of  the  beam  tnat  follows 
the  characteristic  in  the  second  plastic  region  as  shown  in 
Fig.  10.  hue  to  the  slope  of  the  characteristic  the  center  part 
of  this  "cross-section"  has  a coordinate,  denoted  by  s,  that 
differs  from  the  coordinate  x of  the  outer  part  of  the  "cross- 
section".  The  coordinates  are  related  by 

r”.  1 ■ ■ ■ 

s(x)  = x + g(x)  \/l  + 2 sin2f(x)  sin  <p(x) 

= x + \/3  g(x)  sin  f(x). 

The  function  v is  now  determined  at  s by 

t 2 

v[s(x)]  = c - g(x)  \l  1 + 2 sin  f(x)  cos  cp(x) 

= c - g(x)  cos  f(x). 


Hence  the  assumptions  of  continuity  and  of  linear  variation  of 
dx  gives 


d^[s(x) ,z] 


z cos  f(x) 
c - g(x)  cos  f (x7 
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The  equilibrium  equation  is 

3cS  , 

x dx  + xz  _ 0 

3x  ds  6z 


This  determines  t together  with  the  boundary  condition 

xz 

T [s(x),c  - g cos  f(x)]  = X sin  f(x), 

XZ- 


We  obtain 


T 


xz 


[~  sln  f 


+ (1 


(c 


- g cos  f) 


g*cos2f  - cf 1 sin  f i v 

— j 

1 + /3  (g  sln  f ) 


Now  it  remains  to  determine  g(x)  in  terms  of  f(x).  As 
before  we  set  up  the  equation  of  total  moment,  nov/  taken  with 
respect  to  the  point  (s,0).  The  external  moment  is  given  by 
equation  (10)  as 


p s = I — M ( 
F1S  9 „2  pw 


V 


I o 2 M 

-b  PW 


The  plastic  regions  of  the  flanges  contribute  Mn 


Mx  = (a2  - u2)b.jY. 


The  elastic  regions  of  the  flanges  contribute  the  moment  M2 
which  is  made  up  of  two  parts  M2^  and  arising  from  the 

tensile  and  shear  stresses,  respectively 


M2^  = ^ (u  - c)[2u  + c + (u  + 2c)  cos  fjb^Y 
= - — g sin  f[(u  - c)2(l  - cos  f)]!b,Y 


M(2) 

‘2 
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M = f 1 (u  - c)[2u  + c + (u  + 2c)  cos  f] 


-Ig  sinf  [(u-  c)2(l  - cos  f)]‘  l ^Y. 


V3 

The  plastic  regions  of  the  web  contribute  which  also  consists 
of  two  parts,  and  M^2\  from  the  two  stresses,  dx  and  t;xz 

respectively, 

= (2c  - g cos  f)g  cos2f  bY 

= 2gc  sin2f  bY 

M_  = (2gc  - g2cos^f)bY. 

Finally  the  elastic  part  of  the  web  contributes 

Mi,  = ~ (c  - g cos  f)2  cos  f bY. 

f 3 

The  moment  equilibrium  equation  then  becomes 


P„  s = M,  + + M_  + M,  . 

1 12  3 4 


It  is  quadratic  in  g(x).  We  therefore  obtain 


g(x)  = A(x)  - \J A (x)  - B(x) 


where 


2 M 


sin  f + [2u  + 


A(x)  = 


PW 


^3 

“2b  “ 


(19) 


(u  - c)2f’]sin2f 


cos^f 

2 bn 


(20) 


B(x)  = 


2 Plx  ■ Mpf  2 -T 

3c  77 2c  cos  f + ~tr(u  - c)(u  + 2c)  (1  - cos  f) 

M b 

pw 


cos  f 


(21) 
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To  complete  the  analysis  we  must  write  conditions 
ensuring  continuity  of  stresses  at  cross-sections  D and  E.  The 
following  requirements  are  imposed  on  the  values  of  f(x)  and 
f ? (x)  at  D and  E: 


-1 


f(xn)  = cos"  f * (xD)  = 


uDc 


Ur 


Ub  ^ ' °2 

t 

u c 

f(xE)  = cos'1  (-£-)$  f 1 (Xg)  “ 


E 


UE  \f' 


(22) 


(23) 


UE2  ' c2 


The  previous  condition  on  f(x),  (17))  now  takes  the  form 


V3 

2b 


(uF  - c)‘ 


(uD  - 0)' 


Ur 


U 


D 


n*D 


sin  f dx. 


(24) 


E 


Any  function  f(x)  that  fulfills  these  five  boundary  conditions, 
equations  (22),  (23)  and  (24),  and  produces  plastic  regions  of 
the  type  assumed  in  Fig.  7 will  define  a statically  admissible 
system  of  stresses  within  the  framework  of  the  assumptions  we 
have  made.  In  the  numerical  example  one  such  function  f(x)  is 
chosen. 


3.  Beams  of  length  so  that  -L_  > -L  > -i  . 

2 — ~~  0 

The  beams  treated  in  this  part  have  such  length  and 
cross-section  dimensions  that  the  plastic  region  that  started 
from  the  outer  fiber  extends  into  the  web,  and  so  that  the  second 
plastic  region  extends  into  the  zone  FH,  see  Fig,  8.  The  length 
of  the  beam  must  therefore  fulfill 
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l. 


l 


2 o 

— > - > — 

c — c — c 


One  further  condition  should  be  mentioned  namely  the  requirement 
that  the  cross-section  G does  not  fall  near  the  free  end  of  the 
beam. 

The  same  stress  field  as  before  is  used  in  the  zones  BF 
and  GH.  The  "cross-section"  F is  given  by  u^  = a,  i.e., 

(u  - c)2 

(a  - c)(l  - cos  f(x„))  = — — + _2b — j Sin  f dx.  (25) 

UD  b-, 

*JxF 

The  cross-section  G is  determined  by  the  yield  condition  in  the 
web  at  z = c.  It  is 


so 


1 Mrf  2 

“ % (-2£)  . (26) 

M 

pw 

In  the  zone  FG  the  flange  is  entirely  elastic,  l at  the 
web  contains  the  second  plastic  region,  Fig.  8,  The  stresses  in 
the  web  are  determined  as  before  by  the  function  f(x).  The 
expressions  for  the  stresses  in  the  web  remain  unchanged,  so 
that  in  the  plastic  region  of  the  web  along  each  characteristic, 
Fig.  11 

d = Y cos  f,  t;  - JL  sin  f 
x x, 

and  in  the  elastic  region  of  the  web 


G _ 
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dx 


xz 


Z COS  f y 

c - g cos  f 

1 i sin  I + [1 

2 |v/3 


-j  g>  cos2f  - cf 1 sin  f 'j 


2 | 6 ^ I 

2 — ii* 

(c  - g cos  f)  1 + /3  (g  sin  f)  J 


In  the  flange,  which  is  in  the  elastic  state,  we  let 

d have  linear  variation  in  z and  be  continuous  at  the  flange- 
x 

web  junction.  This  gives 

(a  - z)  cos  f + (z  - c)h  _r 


0 = 
X 


a - c 


where  we  denote  g (x,a)  = Yh(x).  The  function  h(x)  will  be  found 
in  terms  of  f(x).  The  equilibrium  equation  (1)  together  with  the 
boundary  condition  t (x,a)  = 0 determines 

XZ 


xz 


= 2^~  :~cT  t(a  + z “ 2c)h»  - (a  - z)f»  sin  f]Y. 


The  boundary  condition  for  t at  the  flange-web  junction  deter- 

XZ 

mines  h(x).  We  have  in  the  flange 


xz 


(x,c)  = ^ (a  - c)(h*  - fT  sin  f)Y  = — Sin  f. 


Hence 


x 

it 


/3  \ 


h = B - cos  f + 


2b  | sin  f dx 
uxG 

V’E  b^(a  - c) 


The  constant  B is  determined  from  the  condition  of  continuity  of 
G at  the  section  G.  We  obtain 

X 


* v.  2b  f sin  f dx 

h = (ft  + 1)  Ji  - - cos  f + 


Pn  C2M 


V 


IY  M 


pw 


V3  bi^a  - °) 


(27) 
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But  the  tensile  stress  is  also  given  at  the  section  F.  This 
imposes  the  following  condition  on  f(x): 


i/3  bx 

-?jj£ — [(a  - c ) ( 1 - cos  f(xp)) 


2 

Pnc  M . 

1 Pfv2n 
J j 


IYM 


pw 


,XF 

j sin  f dx.  (28) 


u'xG 

This  condition  is  equivalent  to  equation  (17)  applying  to  zone 
DE.  It  also  expresses  the  requirement  of  the  equilibrium  of  the 
forces  in  the  x-direction  for  the  segment  of  the  flange  between 


F and  G. 


We  now  determine  the  function  g(x)  by  the  overall  moment 
equilibrium  equation  with  respect  to  the  point  (s,0),  Fig.  11. 

As  before,  the  external  moment  is 


1 


Q 

$ 


2 M 
c p 

7^  m 

pw 


- 1 


The  contributions  to  the  moment  from  stresses  in  the  web  are 
as  given  in  the  previous  section,  namely 

= ( 2gc  - g^cos^f  )bY 

= 2.  (c  - g cos  f)^cos  f bY. 

The  flange  contributes  the  moment  M2,  which  consists  of  two 
parts,  M2^  and  arising  from  the  tensile  stresses  and  the 

shear  stresses,  respectively: 
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i (a  - c)[(a  + 2c)  cos  f + (2a  + c)h]b1Y 

M<2)  = (a  - c)2  (2h  + cos  f)»  bxY 

Mg  = j (a  - c)[(a  + 2c)  cos  f + (2a  + o)h 

+ v^3  (a  - c)g  sin  f(2h  + cos  D’jb-^Y. 

As  before  the  moment  equilibrium  equation 

P-j^s  = Mg  + l' 


determines  g(x)  so  that 


g(x)  = A - \ik2  - B 


(29) 


where  „ - p 2 

c - 3 n/3  F1c 


A(x)  = 


2 M, 


£3L 


sin  f + [2a  + V^.l  (a  - c)2  f«]sin  f 


B(x)  = 


COS^f 

3c2  — i—  - 2c2cos  f - —l(a-c) [(a+2c)cos  f+(2a+c)h] 
M b 

pw  


■(30) 


..  (3D 


COSJf 


As  before,  we  must  check  that  the  stresses  are  continuous 
at  the  cross-sections  F and  G.  The  stresses  at  section  F are 
given  on  both  sides  of  the  cross-section  in  terms  of  f(x). 

Hence  f(x)  must  be  continuous  and  have  a continuous  first  de- 
rivative at  Xp.  Continuity  of  the  shear  stresses  at  F imposes 

t 

one  more  condition,  namely  up  = 0.  In  terms  of  f(x)  this 
condition  takes  the  form 
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M 


ft  (X  ) = - -4  a.  + c JW 


V3 


2 M 


(32) 


At  section  G the  continuity  of  stresses  requires 

2 


f(xQ)  = sin' 


-1  ^ 


2M 


2IY  M 


f (x  ) - - 


pw 


v3c  M 


(33) 


Any  function  f(x)  that  fulfills  the  four  conditions  (28),  (32) 
and  (33)  and  produces  stresses  of  the  type  assumed  in  Fig.  8 
will,  within  our  assumptions,  produce  a statically  admissible 
stress  field.  Note  that  the  position  of  section  F and  the  value 
of  the  function  f(x)  at  F are  not  determined  by  the  stress  con- 
ditions, These  values  may  be  chosen  so  that  the  numerical 
calculations  become  as  simple  as  possible.  It  will  be  convenient 
to  assign  a value  for  x„  at  the  very  beginning  of  the  computa- 

r 

tion.  Thereafter  we  construct  a stress  field  in  the  zone  FG  by 
using  an  f(x)  that  fulfills  the  conditions  (28),  (32)  and  (33). 

We  then  finish  the  problem  in  dealing  with  the  zone  DF,  i.e., 
we  choose  an  f(x)  that  fulfills  conditions  (22),  (24)  and  the 
values  of  f(Xp)  and  f'(Xp)  that  are  obtained  from  zone  FG.  The 
third  numerical  example  in  the  following  section  is  computed 
according  to  this  scheme. 


4.  Numerical  examples. 

In  this  section  an  outline  is  shown  of  the  computations 
that  determine  the  shape  of  the  plastic  regions  for  three  special 
cases  all  applied  to  an  8 WF  40  beam.  The  numerical  values  of 
the  different  quantities  are  given  to  as  many  as  six  significant 
figures,  although  of  course  for  practical  problems  such  precision 
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is  meaningless.  It  is  necessary  in  certain  cases  to  carry  this 
number  of  significant  figures  in  order  to  show  some  of  the  prop- 
erties of  the  analysis. 

We  have  chosen  in  oui  examples  one  of  the  beams  listed 
in  Table  1,  This  table  shows  the  computed  lengths  -t  and 
^2  for  some  beam  sections  which  have  been  subjected  to  labora- 
tory tests  [11] s [12],  [I3]j  and  for  a few  other  common  beam 
sections.  From  these  beams  the  8 WF  40  section  has  been  chosen 
fur  detailed  discussion  because  it  has  a large  value  of  the 
ratios  Mp/Mpv  and  b^/b.  The  cross-section  dimensions  are  given 
in  Table  1.  We  compute  three  examples,  choosing  the  length  of 
the  beam  so  that  one  example  is  provided  for  each  part  of  the 
analysis. 

In  the  first  example  we  consider  the  shortest  possible 
8 WF  40  beam  in  which  the  second  plastic  region  does  not  appear. 
Table  1 gives  -i^/c  = 21,1,  and  we  therefore  take  Z/c  = 21,1  for 
the  first  example.  The  estimate  for  the  limit  load  in  this  case 
is  given  by  equation  (10)  as  P^/M  = 0.9895.  The  beam  contains 
the  zones  BC,  CF  and  FH  where  x^/c  = 20,48,  (Eq,  6),  and 
Xp/c  = 19.02,  (Eq.  8).  The  elastic-plastic  interface  is  computed 
in  the  zones  BC  and  CF  according  to  equations  (5)  and  (7).  The 
results  are  shown  in  Table  2 and  in  Fig.  12. 

In  the  second  example  we  choose  the  length  -i/c  so  that 
Z-^/c  <l/c  < We  take  -t/c  = 12.5356,  (The  actual  choice 

was  to  make  Ug/c  = 0.85  exactly,)  The  estimate  for  the  limit 
load  is  given  by  equation  (10)  as  P1VM  = 0,9715,  The  beam  is 
divided  into  the  zones  BC,  CD,  DE,  EF  and  FH  according  to  Fig.  7. 
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The  elastic-plastic  interface  in  zone  BC  is  given  by  equation 

(5).  The  results  are  shown  in  Table  3»  Equation  (6)  gives 

x /c  = 12. 39b-2,  The  cross-sections  xn/c  and  xp/c  are  determined 
C u hj 

from  un/c  and  uT„/c  in  the  following  way.  Eauation  (12)  deter- 

U Ji 

mines  approximately 

u /c  - 1 = 0.00176;  u,.,/ c - 1 = 0.129. 

D -CJ 

Using  equation  (11)  we  obtain  the  more  accurate  values 
u^/c  - 1 = 0.00175;  u^/c  “ T = 0. lb-2 56. 

Equation  (7)  then  gives 

x / c = 12.3921;  x / c = 11.6367. 

D & 

Cross-section  F is  determined  by  equation  (8)  as 

x _/c  = 11.5074. 

F 

The  zones  CD  and  EF  become  so  narrow  that  the  boundary  values 
of  u/c  are  sufficient  to  determine  the  elastic-plastic  interface 
in  these  zones.  We  therefore  now  focus  our  attention  on  zone  DE, 
In  order  to  construct  a solution  in  zone  DE  we  need  to 
determine  the  function  f(x).  This  is  restricted  only  by  the 
boundary  conditions,  which  are 

f (12. 3921c)  = 0.05917;  f1 (12.3921c)  = - l4.88l/c  (22a) 

f(ll. 6367c)  = 0.50488;  f» (11.6367c)  = - 0.17083/c  (23a) 


and 


All-107 


35 


12.3921c 

'N 

sin  f dx  = 0.3^083c.  (24a) 

11  .'6367c 

We  must  therefore  choose  a certain  form  of  the  function  f(x), 
this  choice  containing  five  arbitrary  constants  so  that  the 
above  conditions  can  be  satisfied.  Note  that  even  if  the  above 
conditions  are  satisfied  the  choice  of  f(x)  may  not  be  satis- 
factory due  to  the  fact  that  the  results  may  not  agree  with  the 

general  type  of  solution  assumed  in  Fig.  7.  E.g.,  the  function 

v(x)  may  become  parallel  at  a point  with  the  z-axis  or  it  may 

become  larger  than  c at  some  point.  Both  cases  must  be  avoided. 
For  the  choice  of  f(x)  in  this  case  we  notice  that  f * (xD)  is 
numerically  large.  This  can  be  taken  account  of  conveniently 
by  choosing  the  f(x)  basically  as  a parabola  with  the  axis 
parallel  to  the  x-axis  and  the  vertex  near  x^.  This  suggests 
the  introduction  of  the  new  variable  y and  also  of  the  square 
root  as  follows: 

y = 12.3921  + k - 2 

c 

and  f (y ) = A + BN/y  + Cy  + Dy2. 

The  five  constants  k,  A,  B,  C and  D are  to  be  determined.  The 
conditions  become  now 

f(k)  = 0.059175  f«(lO  = 14.881 

f (k  + 0.755^)  = 0.504885  f ■ (k  + 0.755^)  = 0.17083 


and 
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k+0 .759+ 

0 

sin  f dy  = 0.34083. 

\J  k 

Prime  denotes  differentiation  with  respect  to  y.  These  con- 
ditions give 

k = 0.0054;  A = - 0,1180;  B = 2.6358$ 

C = - 3.0663;  D - 1.1344. 

The  function  f is  shown  in  Fig.  13.  Equations  (16),  (18)  and 
(19)  determine  u,  9 and  g;  see  Table  4.  The  plastic  regions  are 
shown  in  Fig.  14, 

As  the  third  example  consider  the  shortest  possible 
length  for  v'hich  the  previous  analysis  applies,  i.e., 

Vc  = -co/c  = 10.53^92.  Very  similar  considerations  apply  in 
this  case  as  compared  to  the  foregoing  example.  The  presenta- 
tion is  therefore  shortened.  We  have  the  zones  CD,  DF>  FG  and 

GH  as  in  Fig.  8.  In  this  case  equation  (10)  gives  Pyt/M  ==  0.9605. 

__ 

Equation  (12)  gives  u^/c  - 1 = 0.00119  and  from  equation  (.11) 
u^/c  - 1 = 0,001185.  Equations  (7)  and  (26)  give  x^/ c = 10,53385 
and  xG/c  = 9.29002.  Conditions  (22),  (32)  and  (33)  are: 

f (10. 53385c)  = 0, 0486 5;  f»  (10.53385c)  = - 22.011/c  (22b) 

f'(xF)  = _ 0.33448/c  (32b) 

f (9. 29002c)  = 0,60659;  f«  (9.29002c)  = - 0.l55lVc.  (33b) 


Conditions  (24)  and  (28)  give 
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10.53385c 

sin  f dx  = 0.6835^  (24b)  + (28b) 

9.29002c 

*F 

sin  f dx  = 2,9898  [cos  f(xp)  - 0.77136].  (28b) 

9.29002c 

In  this  example  the  cross-section  F is  arbitrarily  chosen.  We 
take  xp  = 9.5* 

In  zone  FG  we  let 

f = 0.60659  - 0.l55l4y  + Ay2  + By3 

where 

y = x/c  - 9.29002. 

The  conditions  (33b)  are  satisfied  while  A and  B are 
determined  by  (28b)  and  (32b).  The  constants  become  A = 4,1964 
and  B = - 14.679.  The  equations  (18),  (27)  and  (29)  determine 
9,  h and  g in  zone  FG,  see  Table  5. 

In  zone  DF  we  nave  the  following  conditions  on  f(x) 

f (10. 53385c)  = 0.04865;  f‘ (10. 53385c)  = - 22.011/c 

f (9.5c)  = 0.62314;  f«(9.5c)  = - 0.33448/c 

and 

10.53  3 8 5c 

l \ 

sin  f dx  = 0. 56189c. 

u9. 5c 

Now  choose 

f = A + B /y  + Cy  + Dy ^ 


where 
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y = 10.53385  + k - x/c. 

The  above  five  conditions  determine  k = 0.00171?  A = - 0.02903? 

B = 1.93656,  C = - 1.40448  and  D = 0.11020.  Equations  (16), 

(18)  and  (19)  determine  u,  y and  g in  zone  DF?  see  Table  6. 

The  function  f and  the  plastic  regions  for  this  third  example 
are  shown  in  Figs.  15  and  16, 

D.  A second  example  of  an  approximate  statically  admissible 
stress  field. 

The  results  of  the  previous  analysis  in  this  paper  show 
that  for  beams  that  have  a length  ratio  greater  than  -t  /c  the 
reduction  of  the  limit  moment  is  at  most  of  order  10$,  Ex- 
amples of  these  reductions  are  given  in  Table  1 in  the  column 
that  shows  the  value  of  P^£/Mp  for  beams  having  the  length  ratio 
-tQ/c.  The  cases  where  shear  stresses  become  important  therefore 
are  expected  to  occur  for  beams  that  have  length  shorter  than 
£ , A considerable  reduction  in  the  maximum  end  load  of  a 
built-in  I-beam  may  take  place  when  the  beam  becomes  so  short 
that  the  shear  capacity  of  the  web  is  not  sufficient  to  balance 
the  longitudinal  force  in  the  flanges.  A failure?  as  depicted 
in  Fig.  17?  may  then  happen.  As  cross-sections  no  longer  re- 
main plane  and  normal  to  the  center-line  we  can  assume  that  the 
web  carries  a small  amount  of  the  longitudinal  stresses  so  that 
in  the  web  the  shear  stresses  nearly  reach  their  yield-value. 

We  shall  now  construct  a statically  admissible  stress  field  for 
this  case  --  starting  with  the  same  assumptions  as  before  — 
where  this  stress  field  is  based  on  the  failure  of  the  web  in 
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pure  shear.  As  depicted  in  Fig.  18  let  the  stress  field  be 
composed  as  follows: 

Y 


ax  = °s 


d = 


bxY 


xz  >/3 

b(a  - z)Y 


xz 


in  the  web 


in  the  flange. 


v/3  (a  - c)  /3  b1(a  - c) 

The  equilibrium  condition  (1)  is  satisfied  everywhere  and  the 
yield  condition  (2)  is  obviously  satisfied  in  the  web.  At  the 
flange-web  junction  the  shear  stresses  are  in  equilibrium  in  the 
same  sense  as  before,  and  the  longitudinal  stresses  are  discon- 
tinuous. Such  a discontinuity  is  permissible  in  the  construction 
of  a statically  admissible  stress  Cield,  see  [l4],  This  stress 
field  is  valid  when  the  yield  condition  (2)  is  not  violated  in 
the  flanges.  The  yield  condition  in  the  flange  takes  the  form 


bx 


_ 9 


+ 3 


;A'a- z) — l2  < l. 


X I — 


i/3  bx(a  - c)j  |_/3  bx (a  - c)j 

The  left-hand  side  has  its  maximum  value  at  x = <1  and  z = c. 

The  yield  condition  will  therefore  be  satisfied  for  beams  having 
a length  t < ^ where  is  obtained  from  the  following  equa- 

tion : 


b-i- 


VT  b1(a  - c) 


.l2+3^ 


~i2 


V3  bi; 


= 1 


or 


l^/c  = /3  (a/c  - l)\j  (b-j/b)  - 1. 


(34) 


This  length  is  computed  and  shown  in  Table  1,  The  estimate  of 
the  maximum  total  shear  force,  P^,  becomes  in  this  case 
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P = 2bc  •—  + (a  - c)b  ~ = -i-  (a  + c)bY. 
v/3  v'3  v'3 


Therefore 


n M 

rr  ■■■  = ~=  (a/c  +1)  •£./  c —P-v^ 

Mp  % 


(35) 


The  above  stress  field  suggests  a more  general  field  as 
depicted  in  Fig.  19*  We  let 

~ z Y in  the  flange 


dx 

= kl 

£ 

c 

Txz 

n 

H 

dx 

= k2 

X 

c 

y; 

= k2 

We 

assume 

that 

4-  /> 

two 

equilibrium  condition  (1)  is  satisfied  everywhere  in  the  same 
sense  as  before.  The  yield  condition  must  be  checked  in  the 
flange  and  ohe  web.  We  let  1c,  and  k2  be  so  determined  that  the 
yield  condition  is  simultaneously  satisfied  at  the  most  critical 
points  both  in  the  flange  and  in  the  web.  In  the  flange  this 
point  is  x = JL  and  z = c.  The  following  equation  therefore 


determines  k-^: 


(k^c)2  + 3 [kiCa/c  - l)]2  = 1 


k,  = 


1 — • 

\/(<7 c)2  + 3(a/c  - l)2 

In  the  web  the  critical  point  is  x = t and  z = 0.  We  therefore 


choose  k2  so  that 


(k2Vc)2  + 3tk2  + k^o-^/bCa/c  - l)]2  = 1 


*»■/-»  r\ 

~ - + B - A 

kl 


All-107 


41 


where 

A = 3 bx/b  — - aZS-=-L- 

(Vc)  + 3 

('t/C)2  - 3(a/c  - l)2[(b1/b)2  - 1] 

B = i . 

(<e/c)2  + 3 

The  condition  k^  > 0 is  satisfied  but  the  condition  k2  J>  0 shows 
that  B 2 0)  and  this  restricts  the  length  of  the  beam.  This 
restriction  is 

Vc  > /3  (a/c  - 1)  \ (b-j/b)2  - 1 = £ /c.  (36) 

This  stress  field  together  with  the  previous  one  therefore 
applies  to  any  length  of  the  beam.  The  total  shear  force  now 
becomes 

Pp  = k-^(a/c  - l)(a  - c)b^Y  + 2k-^(a/c  - Db^cY  + kpbcY. 

Therefore 

Po't  M j>  f o 

sf-  = nT  I - l)]b1/b  + k2  ] . (37) 

Within  our  assumptions  equations  (35)  and  (37)  provide  a lower 
bound  for  the  limit  load.  The  curve  of  Pp't/Mp  versus  t/c  has 
been  computed  for  the  8 WF  4-0  section.  The  result  is  shown  in 
Table  7 and  in  Fig.  21. 

In  this  stress  field  the  assumption  that  6Z  vanishes  is 
strictly  valid.  This  is  true  because  the  shear  stresses  are 
independent  of  x so  that  the  equilibrium  equation  invol/ing  dz 
is  identically  satisfied.  The  stress  fields  that  we  have  con- 
structed give  a lower  bound  for  the  collapse  load.  In  the  next 
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section  we  shall  obtain  an  upper  bound  which  will  be  applicable 
for  short  beams. 

E.  A kinematically  admissible  velocity  field  for  short  beams. 

We  shall  now  make  use  of  the  second  limit  analysis 
theorem  of  Drucker.  Prager  and  Greenberg  [i4],  which  states 
that  a load  obtained  in  a definite  manner  from  a kinematically 
admissible  velocity  field  is  an  upper  bound  for  the  collapse 
load.  A velocity  field  is  defined  as  kinematically  admissible 
when  it  satisfies  the  imposed  velocity  conditions  at  the  bound- 
ary of  the  body  and  when  the  velocities  vanish  on  those  parts 
of  the  boundary  where  the  surface  tractions  are  not  prescribed. 
In  this  paper  v/e  consider  only  velocity  fields  which  are  con- 
tinuous and  satisfy  the  condition  of  incompressibility 

+ _ o 

9x  9y  9z 

where  u,  v and  w are  the  velocity  components  in  the  x,y  and  z 
directions.  The  upper  bound  P*  of  the  load  is  obtained  by 
equating  the  rate  of  work  produced  by  the  applied  forces  to  the 
rate  of  plastic  work  in  the  body. 

It  seems  intuitively  likely  that  a reasonable  upper 
bound  of  the  collapse  load  should  be  such  that  the  value 
P*i/Mp  = 1 is  exceeded  by  only  a small  amount,  if  at  all.  For 
moderately  long  I-beams  we  have  only  been  able  to  construct 
velocity  fields  that  give  P*t/M  considerably  larger  than  unity. 
We  therefore  confine  attention  to  short  beams  where  a "good" 
velocity  field  is  easily  obtained  and  restrict  ourselves  to  the 
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length  of  the  beams  so  that  P*t/M  < 1.  The  velocity  field 
that  is  chosen  is  depicted  in  Fig.  20.  It  is  symmetric  with 
respect  to  the  middle  surface  OFG.  Part  ABFG  is  rigid  and 
fixed.  The  part  BODE  of  the  flange  is  rigid  and  rotates  with 

the  angular  velocity  w about  the  fixed  point  B*  In  part  ABC 

of  the  flange  let 

u = w = - - x - a + z);  v = 0. 

This  is  a motion  parallel  to  AB,  with  velocities  continuous  on 

AB  and  BC.  In  part  BEOF  of  the  web  let 

u = v = 0;  w = - u>  ( & - x - a + c). 


The  velocity  is  continuous  on  BE  and  the  motion  is  parallel  to 
BF.  As  before  we  use  the  von  Mises  yield  condition  and  there- 
fore also  the  von  Mises  flow  rule.  We  have  constructed  a plane 
velocity  field.  The  rate  of  plastic  work,  W^,  for  such  fields 
is  derived  in  [ 51 j page  214. 

y r • 

w = *—  | rdv 

P 73  ! 

0 

where  dV  is  the  volume  element  and  r is  given  by 


r = 


r/Qu  6w>.2 

L ox  " 3zJ 


+ 


r 9u  , 

+ <3x  J 


We  obtain  both  in  ABC  and  BEOF 


r = 


w. 


The  rate  of  plastic  work  therefore  is 


- ~ Yw[V 


/3 


ABC 


+ VBE0F^ 
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W - M w[b.,/b(a/c  - l)2  + l/o  - a/c  + l]  • 

P /J  pw  1 

The  rate  of  external  work  is  Wg 

W e = P*w(£  - a + c). 


The  limit  analysis  theorem  therefore  gives  the  following  upper 
hound  for  the  limit  moment 


ft  2^/  iVMa/c  - l)2 

^p~  \/3  Mpc  [Vc  - a/c  + 1 

P+'i/M  becomes  unity  for 

/ p 

£/c  = A + y7A  - B 


+ 


where 


and 


A = J [ n/3  Mp/H  - 2b1/b  (a/c  - l)2  + 2 (a/c 

B = *3  M /M  (a/c  - 1). 

^ P pw 


(38) 


1)1 


This  upper  bound,  computed  for  the  8 WF  4-0  section,  is  shown  in 
Table  7 and  plotted  in  Fig.  21  as  a function  of  Z/ c « It  is  seen 
that  the  curve  for  the  upper  bound  load  P*  lies  parallel  to 
and  only  slightly  higher  than  that  for  the  lower  bound  load 
corresponding  to  shear  failure  in  the  web. 


F,  On  the  relations  between  experiments  and  theory. 

The  example  of  the  curve  P^Z/M  versus  Z/c  in  Fig,  21 

C.  p 

clearly  shows  that  the  theory  predicts  a certain  length  ratio 
of  the  beam  below  which  the  load  carrying  capacity  of  the  beam 
decreases  abruptly,  because  of  shear  failure  of  the  web.  This 
critical  length  ratio  is  'i^/c*  In  the  following  comparison  we 


All-107 


*+5 


therefore  look  first  for  the  appearance  of  such  a critical 
length  ratio  in  the  experiments,  and,  when  observed,  we  compare 
the  values  of  the  experimental  and  the  theoretical  length  ratios. 

Secondly  we  want  obviously  to  compare  the  experimental 
and  theoretical  values  of  the  limit  moment  in  each  case.  The 
theoretical  values  are  given  by  equations  (35)  > (37)  (which 
furnish  lower  bounds)  and  (33)  (which  gives  an  upper  bound, 
valid  for  short  beams).  In  the  following  experimental  data  from 
three  sources  are  considered.  In  the  paper  by  Johnston,  Yang 
and  Beedle  [13]  shear  failures  in  the  web  of  I-beams  are  re- 
ported, but  details  concerning  the  length  ratios  and  limit 
moment  values  are  not  given.  Baker  and  Roderick  [11]  made  an 
investigation  of  the  influence  of  shear  on  the  limit  moment  of 
a particular  beam  section  and  presented  the  results  in  the  form 
of  load-deflection  curves  from  which  the  variation  of  limit 
moment  with  length  ratio  can'  be  studied,  Hendry  [12]  has  dis- 
cussed tests  on  a series  of  beams  of  different  cross-sections 
and  gives  the  results  in  a table  where  the  limit  moment  was 
defined  by  deflection  considerations.  Adopting  his  definition 
of  limit  moment  comparsions  are  made  in  this  section  between 
the  experimental  results  and  the  theoretical  predictions  both 
as  to  critical  length  ratio  and  magnitude  of  limit  moment. 

In  the  eighth  Progress  Report  from  Lehigh  University 
[13]  two  examples  of  beams  are  shown  where  the  collapse  of  the 
beams  was  apparently  caused  by  shear  failure  in  the  web.  The 
first  example  is  a 4 I 7,7  beam  that  is  simply  supported  and 
loaded  symmetrically  by  two  concentrated  forces.  Each  end  of 
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the  beam  can  therefore  be  represented  by  the  cantilever  problem 
discussed  in  this  paper.  The  length  ratio  £/c  in  the  test  is 
approximately  2.5.  In  the  theory  equation  (34)  gives  = 2.5, 

as  shown  in  Table  1* 

In  the  second  example  an  8 WF  40  beam  is  loaded  so  that 
moments  numerically  equal  but  opposite  in  sense  act  at  the  ends 
of  the  half -beam,  whose  length  ratio  was  approximately  16* 

Because  of  the  symmetrical  properties  of  the  load  each  quarter 
of  the  span  can  be  represented  by  the  cantilever  beamB  The 
effective  length  ratio  therefore  is  -t/c  =8,  For  this  case  we 
compute  £^/c  = 6.0. 

The  ratio  between  the  collapse  moment  and  the  fully 
plastic  moment  is  not  stated  for  these  two  examples.  The  tests 
show  that  shear  failure  in  the  web  is  important  for  beams  hav- 
ing the  above  cross-sections  and  the  effective  length  ratio 
approximately 

Simply  supported  beams  acted  upon  by  two  concentrated 
loads,  symmetrically  spaced,  have  been  investigated  and  re- 
ported by  Baker  and  Roderick  [ 11] « The  beam  is  a British  H- 
section  with  dimensions  I-4-11  x If".  The  distance  between  the 
loads  was  kept  constant  and  the  span  was  varied.  Different 
effective  length  ratios  of  the  corresponding  cantilever  beam 
problem  were  therefore  obtained,  varying  from  9*4  (Test  No,  SR 
21)  to  2.2  (Test  No.  SR  11),  The  value  of  £.^/c  is  6,7.  The 
load-deflection  curves  are  given.  For  length  ratios  greater 
than  3.9  the  curves  show  a well  defined  "bend"  where  the  slope 
of  the  curve  rapidly  decreases.  The  load  at  these  "bends"  is 
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taken  as  the  collapse  load  in  Table  8 and  Fig.  22.  The  table 
and  the  figure  show  that  in  the  experiment  when  the  length  ratio 
is  smaller  than  5*5  the  collapse  moment  decreases  rapidly  with 
further  decrease  in  length  ratio.  In  the  theory  we  have 
'L/c  = 6,7  below  which  decreases  rapidly. 

In  the  paper  by  Hendry  [12]  the  influence  of  shear 
forces  in  simply  supported  beams  of  various  I-sections  is  dis- 
cussed. In  these  tests  a concentrated  load  acts  on  the  mid- 
point of  the  beam.  A series  of  tests  is  also  included  in  which 
the  influence  of  shear  stresses  in  the  beam  of  a portal  is 
investigated.  The  value  of  M/Mp  is  tabulated  in  the  paper, 
where  the  moment  M is  obtained  from  the  load  at  which  the 
deflection  of  the  beam  equals  1/50  of  the  span.  These  values 
are  shown  in  Table  9 and  Figs,  23  and  24  together  with  the 
theoretical  values  of  P't/Mp,  Except  for  the  portals  and  the 
beam  3b  we  sea  that  M/M  is  in  the  neighborhood  of  and  is  great- 

ir 

er  than  W e notice  moreover  that  for  length  ratios 

shorter  than  -03/°  the  value  M/Mp  is  considerably  less  than 
unity. 

The  portals  B1  - B5  are  made  up  of  a 3”  x 1"  beam  and 
of  1-4-"  x 1"  stanchions.  The  ratio  of  the  fully  plastic  moment 
of  the  beam  to  that  of  the  stanchion  is  approximately  7»3  to  1, 
The  effective  length  ratio  Z/c  is  therefore  taken  to  be 
—j  7*3.  — (Z/c)*  - 0,88  (Z/c)*  where  (Z/c)*  is  the  length  ratio 
between  the  load  and  the  end  of  the  beam.  The  results  in  Table 
9 show  that  the  portals  did.  not  develop  the  moment  P2Vm  . 
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In  [12]  the  large  reduction  of  the  limit  moment  is  explained 
by  the  influence  of  the  axial  stresses  in  the  columns. 

In  general  we  see  that  the  theory  and  experimental 
tests  agree  rather  well,  with  very  few  exceptions.  For  this 
comparison  between  the  experiments  and  the  theory  we  have  used 
the  results  obtained  by  the  second  type  of  stress  field.  This 
is  done  because  the  length  ratio  is  in  these  cases  so  short  that 
the  first  type  of  stress  field  does  not  apply. 

G.  Summary  and  conclusions. 

The  two  types  of  stress  field  in  this  paper  each  give 
an  approximate  lower  bound  for  the  true  collapse  load  of  the 
plastic  bending  of  a cantilever  I-beam.  The  approximations 
consist  in  both  cases  of  a two-dimensional  formulation  of  the 
problem  while  in  the  first  type  of  stress  field  the  further 
assumption  is  made  that  the  a stress  Is  negligible.  The  veloc- 
Ity  field,  which  is  applied  for  short  beams,  gives  results 
which  are  upper  bounds  to  the  true  collapse  load.  The  velocity 
field  satisfies  all  the  requirements  of  the  appropriate  theorem 
of  limit  analysis.  The  upper  and  lower  bounds  differ  by  small 
amounts  for  short  beams.  For  long  beams  the  lower  bounds  fur- 
nished by  both  types  of  stress  field  tend  toward  the  result 
Pi/Mp  = 1. 

The  ratio  P2't/Mp , where  is  "the  lower  bound  on  the 
collapse  load  defined  In  Section  D,  can  be  used  to  estimate 
the  limit  moment  M0  = in  the  relation  to  the  fully  plasitc 
moment  which  is  appropriate  for  pure  bending  (or  very  long 
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cantilevers).  We  note  that  the  curve  of  as  °f 

l/c  has  two  distinct  parts,  as  exemplified  by  the  curves  of 

Figs.  21  - 24.  For  l/c  < l^/c  it  is  a straight  line  through 

the  origin  with  the  slope  a/ 9.  ± . For  l/c  > l^/c  the  value 

Mp/Mpv  3 

of  P^-t/Mp  increases  asymptotically  to  the  value  P2^/^p  = 1» 

The  reduction  of  the  limit  moment  at  the  length  ratio  l^/c  ma^ 
be  as  high  as  30%>  for  the  beams  mentioned  in  Table  1.  The 
curves  show  that  a much  larger  reduction  of  the  limit  moment 
takes  place  when  the  length  ratio  is  less  than  -L^/c,  but  for 
beams  having  a length  ratio  greater  than  -t^/c  the  fully  plastie 
moment  is  generally  reduced  by  a relatively  small  amount.  Table 
1 also  shows  that  the  value  of  l^/Q  in  some  cases  is  so  short 
that  in  the  practical  problem  the  effects  of  the  details  of  the 
end  constraint  or  of  the  loading  are  very  important.  It  is 
therefore  suggested  that  beams  should  not  be  used  whose  length 
ratio  is  smaller  than  1^/°  or  which  are  so  short  that  the  end 
effects  are  important.  For  example,  if  l/c  < 6 treatment  as  a 
"beam"  in  the  conventional  sense  is  hardly  appropriate. 

The  influence  of  buckling,  particularly  of  the  web  or 
the  flange,  (or  both),  is  disregarded  in  this  paper.  This  may 
be  important  for  certain  short  beams,  as  indicated  in  [13], 
and  should  be  further  investigated. 
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Table  3j»  The  function  u/c  in  zone  EC  for 
8 WF  40  section; -i/c  = 12.5356 


x/c 

V 

» 

f 

r~ 

u/c 

9 

j 

; g/c 

t 

/"  "2“ 

g/c  v J 1+2 sin  f 

12,3921 

0.0054 

0.05917 

1.00175 

5.9° 

jo. 0000 

0.0000 

12.3 875  ' 

0.01 

0.11504 

1.0036 

11.3° 

jo. 0011 

0.0011 

12.3775 

0.02 

0.19390 

1.0056 

I8.80 

j 0.0079 

I 

0.0082 

12.3675 

0.03 

0.24759 

1.0072 

23.6° 

1 

' 0.0166 

0.0176 

12.3575  ■ 

0.04 

0.28835 

1.0087 

27.2° 

*0.0250 

0.0270 

12. 3475  ! 

0.05 

0.32094 

1.0101 

29.9° 

0 

• 

0 

UJ 

-r 

-r 

0.0377 

12.3275  : 

0.07 

0.37032 

1.0128 

33.9° 

; 0.0497 

0.0573 

12.2975 

0.1 

0.42025 

1.0169 

38.5° 

0.0666 

0.0769 

12.2475  | 

l 

0.15 

0.46845 

1.0241 

41.20 

0.0814 

0.0966 

”7 

12.1975  ! 

0.2 

0.49291 

1.0319 

42.9° 

i 0.0852 

0.1025 

12.0975 

0.3 

0.50793 

1.0500 

44.0° 

| 0.0744 

0.0903 

11.9975 

0.4 

0 . 50408 

1.0711 

43.7° 

0.0536 

. 

0.0649 

11.8975 

0.9 

1 0.49631 

1.0940 

43.2° 

; 

: 0.0317 

0.0383 

11.7975  ; 

0.6 

0.49235 
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42.9° 

' 0.0139 

0.016  7 

11.6975 

] 

0.7 

; 0.49679 
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; 0.0031 

O.OO38 
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0.50488 
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1 
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Table  7.  The  values  of  P,-t/M  , 
1 P 

P^t/Mp  and  for  8 WF  40  section 


g/c  jg/c  \Jl  + 2sin^f 


10  . 53335  0.00171 
10.52556  0.01 


0.02 


0.03 


10. 4-9556  0.04 
10.48556  0.05 


10.46556  0.07 
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0.21675  jl. OO527! 20.9° ;0. 0086 


. I i 
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Table  6. 

The  functions  f,  u/c.  © 

and  e/c 

in  zone  DF 

for  8 WF  40  section;  Uz  = 10.53492 
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Test  No, 

i l 

i 

i 

I l/c 

i 

r.""—  — — ■■■• 

■ 2W 

| 

i 

! VJt/c  j 

l 

i P*l/H 

P 

SR  21 

■ 4.25 

i 

i 

9.44 

| 2.05 

1 

9.68  i 

0.918 

> 1 

20 

j 3.75 

! 8.33 

2.30 

9.58 

0.904 

l 

i 

> 1 

19 

! 3.25 

7.22 

1 2.65 

i 9.57 

0.881 

0.949 

f 

22 

3.00 

6.67  | 

2,85 

i 9.50 

0.857 

| 0.891 
j 

17 

; 2.75 

6.11  ! 

i 

3.0  5 

9.32  ; 

0.785 

! 0.831 

16 

2.50 

1 5.56  | 

3.30 

9.17  1 
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| 0.773 

15 

! 2.25 

■ 5. CO 

3.15 

| 7.88  ' 

0.643 

0.715 

10 

2.00 

! 4.44 

3.15 

6.99 

0.571 

0.656 

14 

1.75 

i 

■■i  - . , 

| 3.&9 

: i 

; 

3.30 

! 

:■  6.42 

i 

i 

1 

0.500 

i 

1 

| 0.600 

Table  8.  The  tests  by  Baker  and  Roderick 
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1 

Beam  No.  j 

i 

Section 

7^ 

i 

l/c  1 

I 

p2Vmp 

p*Vm 

p 

M/Mp 

la 

i 

1 

15 

9.1  i 

0.973 

> 1 

0.95 

b 

4-"  x 3"  R.S.J. 

12 

7.3  ' 
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> 1 

1:00 

c 
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9 

I 
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> 1 

1.00 

I 

d 

6 

3.6  ■ 
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0.83 

e 

' 
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0.5H 
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2a 

i 

12 

77 
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b 
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c 
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5 
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d 

1 

3 
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9 
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b 
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4 
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f 
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i 

10 

i 

* 
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i 

Table  9.  The  tests  by  Hendry 


. 15  and  16.  The  function  fix)  and  the  plastic  regions  for 


Fig.  19.  The  second  type  of  stress  field 
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Fig.  21.  The  theoretical  values  of  the  limit  load  for  8WF40sectloh 
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